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Abstract 

If the recently-discovered charmonium state X(3870) is a loosely-bound molecule of the charm 
mesons D° and D*° or D° and D*°, it can be produced in e + e~ annihilation at the T(45) resonance 
by the coalescence of charm mesons produced in the decays of B + and B~ mesons. Remarkably, in 
the case of 2-body decays of the B mesons, the leading contribution to the coalescence probability 
depends only on hadron masses and on the width and branching fractions of the B meson. As the 

1/2 

binding energy Ej, of the molecule goes to zero, the coalescence probability scales as E b \og{E^). 
Unfortunately, the coalescence probability is also suppressed by two powers of the ratio of the 
width to the mass of the B meson, and is therefore many orders of magnitude too small to be 
observed in current experiments at the B factories. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40. Gx 
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The recent unexpected discovery of a narrow charmonium resonance near 3.87 GeV 
presents a challenge to our understanding of heavy quarkonium. The new charmonium state 
X(3870) was discovered by the Belle collaboration in electron-positron collisions through 
the 5-meson decay B ± — > K ± X followed by the decay X — > J/ipir + iT~ [1]. Its mass was 
measured to be Mx = 3872.0 ± 0.6 ± 0.5 MeV [1]. It is narrow compared to other char- 
monium states above the threshold for decay into DD: the upper bound on the width is 
T x < 2.3 MeV. The discovery has been confirmed by the CDF collaboration who observed 
X through J/ipn + n^ events in proton-antiproton collisions and measured its mass to be 
M x = 3871.4 ± 0.7 ± 0.4 MeV [2]. 

The proposed interpretations of the X(3870) include a D-wave charmonium state with 
quantum numbers J PC = 2 or 2~ + , an excited P-wave charmonium state with J PC = 1 ++ 
or 1 H , a "hybrid charmonium" state in which a gluonic mode has been excited, and a D°D*° 
or D°D*° molecule [3-12]. This last possibility is motivated by the fact that the X(3870) is 
extremely close to the threshold 3871.2 ± 0.7 MeV for decay into the charmed mesons D° 
and D*°. 

If the X(3870) is an S-wave D D*°/D°D* molecule, its binding energy is smaller than 
any other hadron that can be interpreted as a 2-body bound state of hadrons. For two 
hadrons whose low-energy interactions are mediated by pion exchange, the natural low- 
energy scale for the binding energy of a molecule is m^/(2m), where m is the reduced mass 
of the two hadrons. The natural low-energy scale for a D°D*° molecule is about 10 MeV. The 
measurements of the mass of the X(3870) imply that its binding energy (which is positive 
by definition) is E b = -0.5 ± 0.9 MeV. Thus E b is likely to be less than 0.4 MeV, which is 
much smaller than the natural low-energy scale. 

The tiny binding energy of the X(3870) implies that the D°D*° scattering length a is 
unnaturally large compared to the natural scale l/m n . The molecule therefore has properties 
that depend on a but are insensitive to other details of the interactions of D° and D*°, 
a phenomenon called "low-energy universality." The universal prediction for the binding 
energy is 

E b = m D + m D » - M x ^ ? , (1) 

2m DD *a 2 

where m^D* = m D om D *o / {m D o + m D *o) is the reduced mass of the two constituents. The 
universal prediction for the momentum space wavefunction of the D°D*° or D°D*° is 

(Svr/a) 1 / 2 

^ ~ W+lp \k\«m n , (2) 

where the normalization is J d 3 k/(27r) 3 \ifj(k)\ 2 = 1. This wavefunction has been exploited 
by Voloshin to calculate the momentum distributions for the decays X — > D D ti and 
X — > D°D°'y [6] . The universal prediction for the amplitude for elastic scattering of D° and 
D*° in the center-of-momentum frame with total energy E is 

A[D°D*° — > D°D*°] ~ 8 ; m " m r„ h \k\«m n , (3) 

m DD *(-l/a-i\k\) 

where \k\ = [2m£,£)*(E — vtld — m/).)] 1 / 2 . Other consequences of low-energy universality 
have been discussed in Ref . [9] . One consequence is that as the scattering length a increases, 
the probabilities for components of the wavefunction other than D°D*° or D°D*° decrease 
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as 1/a. In the limit a — > oo, it becomes a pure molecular state. If it has charge conjugation 
C = ±, the state is a superposition (\D D*°) ± \D°D*°))/V2. 

One of the challenges for the interpretations of the X(3870) as a D D*°/D°D* molecule 
is to understand its production rate. The large D°D*° scattering length implies that a D° 
and D*° with relative momentum small compared to have very strong interactions. One 
way to produce X is therefore to produce D° and D*° with small relative momentum which 
then undergo a final-state interaction that binds them to form X. An example of such a 
process is the decay of T(AS) into B + B~, followed by decays of the B + and B~ into states 
containing D° and D*°, respectively. There is a small probability that the D° and D*° will 
be produced with sufficiently small relative momentum for them to coalesce into X. In 
this paper, we calculate the leading contribution to the coalescence probability in the case 
of 2-body decays of the B + and B~. We show that the coalescence probability scales as 

1 /2 1 /2 

E b \ogE b as the binding energy of X goes to 0. Remarkably, the coefficient of E b logE b 
depends only on hadron masses and on the width and branching fractions of the B meson. 

We consider the decay T(4S f ) — > Xhh', where h and h! are light hadrons. This process 
can proceed via the decay T(4S f ) — > B + B~, followed by the 2-body decays B + — > D°h 
and B~ — > D*°h', followed by the coalescence D°D*° — > X. This process can also proceed 
through a second pathway consisting of the 2-body decays B + — > D*°/i and B~ — > £>°/i' 
followed by the coalescence D*°D° — > X. In principle, these two pathways can interfere. 
However, we shall see that the momentum configurations of X, h and h' are completely 
determined by the masses of the hadrons and they are different for the two pathways. Thus 
there is no interference between the amplitudes. 

The decay process T(AS) — > Xhh' is complicated because there are many relevant energy 
and momentum scales and they range over many orders of magnitude. The mass My of the 
T(4S) is larger than the binding energy E b of X by more than 4 orders of magnitude and 
the width of the B meson is smaller by about 10 orders of magnitude. We expect the 
rate for this decay to go to in the limit E b — > (with E b 3> r B ), because the D° and 
D*° must have relative momentum k of order {mE>D*E b ) 1 / 2 m 1/a in order to bind to form 
X and such small relative momentum accounts for a decreasing fraction of the total phase 
space available to the D° and D*°. Our calculation shows that there are contributions to 
the rate that scale as E\J 2 . They include contributions from relative momentum k ranging 
from the scale 1/a to the scale m n . The contributions from k ~ 1/a are constrained by low- 
energy universality, and we expect these to be calculable in terms of the scattering length. 
The contributions from k ~ m n necessarily involve the full complications of low-energy 
QCD. Fortunately we find that there is a logarithmic contribution coming from the range 
1/a C C m n , which is also governed by low-energy universality. This logarithmic term 
dominates in the limit E b — > 0. The logarithmic term in the decay rate for T(4S) — + Xhh' 
is calculated in Appendix A. The branching ratio that measures the coalescence probability 
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for D°D*° -> X or D*°D -> X is 

r[T(45) -> X/i/i'] 



T[T(45) -> D°D*°hh'} + r[T(4S) -> D*°D°hh'} 



m DD *M r (M^ — 4m 2 B y/ 2 \m DD * J \2m DD *E b J \m B 

^ \ 1 ' 2 {m B ,m D ,m h ) \ L ' 2 (m B ,m D *,m h >) 

x -> Z)°/i]i3[i? _ -> D* //]) ~* , (4) 

where J(Mt, • • •) is the function of hadron masses given in (A21) and the function 
A(mi, 7712,7713) is given after (A4). The sum in the numerator and the denominator is over 
two terms, the one shown and a second term obtained by replacing D° and D*° by D*° 
and D°. Notice that the expression (4) depends only on hadron masses and on the width 
T B and branching fractions of the 5-meson. If h and h! are each others antiparticles such 
as 7r + and 7r~, the branching fractions cancel between the numerator and denominator. If 
we take the binding energy of X to be E b — 0.1 MeV, the branching ratios in (4) for the 
cases hh' = (tt+tt- , p + p~ , K + K~ , K* + K*-) are (1.1,1.3,1.2,1.3) x 1(T 24 . For any other 
combination of h = (n + , p + , K + , K* + ) and h' = (ir~ , p~, K~ , K*~), the branching ratio in 
(4) depends on B + branching fractions but it is in the range from 1.2 to 1.4 x 10~ 24 . 

We can get a simple expression that can be used to estimate the order of magnitude 
of the branching ratio by neglecting the light hadron masses m h and m' h , and making the 
approximations m D * — m D <C M x and M x <C m B . In this limit, the function J(pq,p'q) 
given by (A21) approaches 



T( ' \ lxMx 
J{PQ,Pq) — m | (M 2 _ 4m |)i/2- 



(5) 



The branching ratio in (4) then reduces to 

T[T(AS) -> Xhh'} 
T[T{AS) -»■ D°D*°hh'} + T[T{AS) -> D*°D°hh' 



ln 2 m 2 B M x ( 8E b \ 1/2 ( 2m 2 \ ( T B 



log TTET -T- • (6) 



M T (M 2 - 4m|) V^x/ \M x E b J \m B 

If we take the binding energy to be 0.1 MeV, this estimate for the branching ratio is 6.3 x 
10~ 25 , which is within a factor of 2 of the more accurate results calculated using (4). This 
estimate applies equally well if h or h! is replaced by a multiparticle state of light hadrons or 
a lepton pair whose invariant mass is small compared to m B . We conclude that the inclusive 
branching fraction for T(AS) — > X(3870) via this coalescence mechanism is about 24 orders 
of magnitude smaller than the product of the inclusive branching fractions for B + — > D° 
and B+ -> D*°. 

We have calculated the leading contribution to the probability for charm mesons produced 
by the decay of T(45') to coalesce into X(3870). Remarkably, this coalescence probability can 
be expressed completely in terms of hadron masses and the width and branching fractions 
of the B meson. Unfortunately there is a suppression factor of (T B /m B ) 2 that makes the 
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B+(p + £l 

T ( p ) O LJ^X(Q) 



■h'{k') 

FIG. 1: Feynman diagram for the amplitude of T(4S) — > Xhh' via the first pathway. 

rate for T(4S') — > Xhh' many orders of magnitude too small to be observed at the current 
5 factories. 

This research was supported in part by the Department of Energy under grant DE-FG02- 
91-ER4069. 

APPENDIX A: CALCULATION OF THE RATE FOR T(45) > Xhh' 

In this appendix, we calculate the rate for the decay T(4S) — > X(3870) + h + h', 
where h and hi are light hadrons. This decay proceeds through two pathways: the de- 
cay T(4S) -> B + B- followed by the 2-body decays B + -> D°h and B~ -> D* // followed 
by the coalescence process D°D*° — > X, and the pathway obtained by replacing D° and D*° 
by D*° and As we shall see, the two pathways do not interfere. The amplitude for the 
first pathway can be represented by the 1-loop Feynman diagram with meson lines shown 
in Fig. 1. We denote the T(45) simply by T. The momenta of T, X, h and h! are P, Q, 
k and k', respectively. The momenta of the virtual B + , B~ , D° and D*° mesons are p + £, 
p' — I, q + £ and q' — £, respectively, where I is the loop momentum. Momentum conservation 
requires P = p + p' , p = k + q, p' = k' + q' , and q + q' = Q. We can choose q and q' to be 

<t" = (A1) 

so that they are on the mass shells of the D and D*°: q 2 = m 2 D and q' 2 = m 2 D ,. Since 
the binding energy of X is so tiny, the momenta (Al) are consistent with the constraint 
q + q' = Q. 

The decay rate can be written as 

d 3 k d 3 k' d 3 Q 



T[T -> Xhh'} = -4- [ 
1 J 2M T J 



(27r) 3 2fc (2tt) 3 2^ (2vr) 3 2g 

x (2n) 4 5 {4 \P-k-k' -Q)\A[T ^ Xhh']\\ (A2) 
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The amplitude for the decay through the first pathway is 



Ai[T -> Xhti] = J -0^ A[T -> £+£-].A[£ + -> D°h]A[B- -> 2)* //]„4.[2) 2)* -> X] 



x 



(p + f) 2 -m B + im B T B {p' - I) 2 - m 2 B + im B Y B 



X (q + £f -m 2 D + ie (q> - £f - m 2 D * + ie (A3) 

The rate depends crucially on the width T B ot the B meson, so the effect of the width must 
be included in the propagators of the B + and B~ . 

Since the loop integral is dominated by very small momenta, we can neglect any momen- 
tum dependence of the amplitudes A for T — > B + B~, B + — > D°h and B~ — > 2}*°/i'. They 
can be approximated by the amplitudes for the on-shell decays. For example, the amplitude 
for B + — > D°h is determined by the branching fraction for that decay: 

B[B+ - 2)°/i] = J- \A[B+ - D Q h}\ 2 ^f^ , (A4) 
L J 16tt 1 L J| m|r B v 1 

where X(x, y, z) = x 4 + y 4 + z 4 — 2(x 2 y 2 + y 2 z 2 + z 2 x 2 ). The amplitude for the coalescence 
process D°D*° — > X can be deduced from the amplitude (3) for the scattering process 
D°D*° — > D°D*°. This amplitude has a pole in the total energy 22 at the mass M x = 
m D + m D * — E b , where E b is the binding energy given by (1). Its behavior near the pole is 



m 



2 

DD* 



E - (m D + m D * - E b ) ' 



The residue is proportional to the square of the amplitude for D°D*° —> X: 

A[D°D*° - X] = ( ^Mxm D m D .\ 1/2 ^ 
V m DD ,a J 

Our strategy is to manipulate the decay rate into a form that includes as a factor the 
differential decay rate for T — > B + B~ . The first step is to integrate over the component £ 
of the loop momentum. The dominant contribution to the integral over £ in (A3) comes 
from the particle poles in the propagators for the 5° and D*° mesons: 

d J± 1 1 = i ( A7) 

2tt (q + if -m 2 D + ie (q' - t) 2 - m 2 D * + ie AE D E D * {E D + E D * - Q ) ' V ; 

In the rest frame of X, the energies are Ed = (m D + I 2 ) 1 / 2 , Ed* = (m 2 D * + I 2 ) 1 / 2 and 
Qo = M X - Expanding the denominator to lowest order in £ and 1/a, (A7) becomes 

1 _ moo* 1 (Ag) 



AE D E D * {E d + E D * - Qo) 2m D m D * I 2 + 1/a 2 ' 
which is proportional to the momentum-space wavefunction ip(£) in (2). 
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If not for the loop momenta, the product of the B + propagator in (A3) and its complex 
conjugate could be expressed as a Breit-Wigner resonance factor. If we take into account 
the loop momenta, that product can be approximated as 



(p + £) 2 — m? B + iniETs \ (p + (- 1 ) 2 — rn 2 B + irriB^B 
1 im B r B 



p ■(£-£') + {£ 2 - £' 2 )/2 + im B T B (p 2 - m 2 B ) 2 + (m B r B ) 2 

^ k.{Z-Z>)+im B T B ™tf-m\)- (A9) 

In the second line, we expressed the product of propagators in terms of a difference between 
propagators and took the limit £ — > and £' — > in that difference to get a Breit-Wigner 
resonance factor. In the third line, we took the limit Tb — > in the resonance factor to get 
a delta function. We also used the relations £ 2 = —2q ■ £ and £' 2 = —2q ■ £', which follow 
from the fact that q, q + £ and q + £' are all on the mass-shell of the D° meson. The inner 
product k • (£ — £') in the denominator of (A9) can be expanded in powers of the momenta £ 
and £'. The terms q ■ £ and q ■ £' are already second order and could be neglected. However, 
for the purpose of evaluating the integral over £, it is more convenient to use the fact the 
q ■ £ and q ■ £' are second order to replace k^ by a vector kg whose /i — component vanishes 
in the rest frame of X: 

fcg = F - ^-Q». (AlO) 

The expression for q in (AlO) implies kg = Pg. Thus our approximation for the inner 
product in the denominator of (A9) can be written 

k.(£-£')~p Q . {£-£'). (All) 

We can integrate in the momentum p of the B + using the identity 

The integral over p 2 can be evaluated using the delta function in (A9). After similar ma- 
nipulations involving the B~ momentum, our decay rate through the first pathway can be 
written 



I\[T -> Xhti] = J dY[T -> B+B~] \A[B + -> D°h]\ 2 \A[B~ -> D*°ti] 

nMx - fl d3Q (MS) 



m D m D * 



where dT[T — > B + B ] is the differential decay rate for T into B + and B with momenta p 
and p'. In the rest frame of X, the factor / in (A13) is given by the integral 

f dH dH' l 1 l l__ 

J (2tt) 3 (2tt) 3 -pg ■ (£ - £') +iep Q -(£- £') + ie £ 2 + 1/a 2 £' 2 + l/a 2 ' { } 
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We have replaced the terms iniB^B in the propagators by ie, because the integral has a 
well-behaved limit as Tb — > 0. 

In order to evaluate /, we first combine both the denominators £ 2 + m 2 and £' 2 + m 2 where 
m = 1/a and the denominators —pQ-(£—£')+ie and p'g- (£—£') +ie into squared denominators 
using Feynman parameters x and z. We then combine the squared denominators using an 
integral over y: 



-p Q ■ (£ - £') + ie p' Q ■ (£ - £') + ie £ 2 + m 2 £> 2 + m 2 



i x ,1 
dx -^—r - — — — I dz 



o 



[x£ 2 + (1 - x)£' 2 + m 2 ] 2 J [C(z) ■ {£ - £') + ie} 2 

= Q l dX l dZ Jo ^ [x£ 2 + (1 - x)£' 2 + m 2 + yC(z) • {£ - £') + ief ' (M5) 

where C(z) = — zp Q + (1 — z)p' Q . The integrals over the momenta £ and £' can be simplified 
by first shifting them and then rescaling them by factors of x and 1 — x: 



I 



dH dH' 



(2tt) 3 (2tt) 3 [x£ 2 + (1 - x)£' 2 + m 2 + yC ■ (£ - £') + ie] 4 
d 3 £ dH' 1 



/ 



= x 



(2tt) 3 (2vr) 3 [x£ 2 + (1 - x)£' 2 + m 2 - y 2 C 2 /(4:r(l - x)) + ie] 

i/2/i _ T r 3 / 2 / ^ (A]6) 

1 J J {2^f{2nf\£ 2 + £' 2 + m 2 -y 2 C 2 /{Ax{l-x))+ieY' K ' 



The integrals over y and then x can now be evaluated analytically. The resulting expression 
for I is 

f 1 A i r i 

" 47r 7 ^C(^) 2 J (2vr) 3 (2vr) 3 [£ 2 + £' 2 + 1/a 2 ] 3 ' (Ai7j 

The integral over £ and £' has a logarithmic ultraviolet divergence. It can be regularized by 
subtracting the integral with m — 1/a replaced by an ultraviolet cutoff A m: 

d 3 Z d 3 £' 1 , 1 , A /A1oN 

(™-> A ) = ^kg-. (A18) 



(27r) 3 (27r) 3 [£ 2 + £' 2 + m 2 ] 3 v ' 64vr 3 to m' 

The appropriate choice for the cutoff is A = m n , since the region of validity of the expression 
(2) for the wavefunction of X is \k\ <C m n . The integral over z in (A17) can be expressed 
in a manifestly Lorentz-invariant form. Our final expression for the integral / in (A14) is 

1 = i6"^ log(/im ^ J (PQ>p'q)' (A19) 
where J(pq,Pq) is a Lorentz-invariant function of pq and p'q defined by the integral 

J(PQ,P' Q )= [ dz- ~ l (A20) 

Jo [zp Q - (1 - z)p Q \ 
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The integral can be evaluated analytically: 

JM = D^j + arCtan ^(S) ' <A21) 

1 /2 

where pq, p' Q and Pq are all defined by (A10) and D(pq,p'q) = [pqPq — (pq ■ Pq) 2 ] ■ In 
the rest frame of X where pq and p' Q are spacelike, D(pq,p'q) = \pg x p'q\. 

It remains only to evaluate the integrals over k, k! and Q in (A13). The integral over k 
can be evaluated by using the identity 

/rl 3 k 
^-^-(27r) 4 5( 4 )(p -q-k) = 2t:5((p - q) 2 - m\\ (A22) 

and similarly for k! . The two remaining delta functions can then be used to evaluate the 
integral over Q: 

d 3 Q n_r//__ _.\2 „_2\ n_ £•//__/ J\2 ...I2\ 



I j2^Q 27r6{iP ~ qf - ml) 27rW - q ' )2 ~ ^ 



4m D m D *M r (M 2 - 4m 2 B y/ 2 ' 

(A23) 

Our final expression for the decay rate through the first pathway is 

r\[T -> xhti] =r[T -> p+p-] -> z) /i]| 2 |„4[p- -> p*°/i']| 2 

M xJ(Pq,Pq) \og(am w ) 
Q47rm DD *M r (M 2 -4m 2 B y/ 2 am DD * ' 1 J 

The function J{pq,Pq) given explicitly in (A21) is a Lorentz-invariant scalar function of the 
momenta p, p' and Q. The Lorentz invariants are p 2 = p' 2 = m 2 B , Q 2 = M\ and 

r\ M * (2,2 2\ 

V • Q = 7^ \ m B + ™ D - ™ h ) , 

' r\ Mx (2,2 2\ 
P - Q = 7. \ m B + m D* - m h> ) , 

p-p'=^(M 2 -2m 2 B ). 

The decay rate through the second pathway is obtained by replacing D° and D*° by D*° 
and D°. The inner products between the momenta k, k! and Q of the final-state particles 
are completely determined by the hadron masses. Since these inner products are different 
for the two pathways, they produce distinct momentum configurations and they therefore 
cannot interfere. 

The rate for the decay T -> P+P~ followed by the decays B + -> D°h and B~ -> D*°h' 

is 

ri[T -> D°D*°hti] = L[T -> P + P~] P[P + -> P°/i] tf[P~ -> P*°/i']. (A25) 

The expression for the branching fraction P[P + — > D°h] is given in (A4). The amplitude 
A[B+ -> D°/i] in (A24) can be eliminated in favor of P[P+ -> D°h]. 
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